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Abstract. Using an experimental based optimization approach, this study investigated whether two specific 
human movements, i.e. pointing tasks and squat-jumps, could be controlled with a reduced set of kinematic 
parameters. Three models were proposed to model the evolution of joint angles as sigmoid shaped curves. 
304 pointing tasks and 120 squat-jumps performed by 22 and 13 subjects respectively were used to quantify 
the accuracy of the models. The parameters of the sigmoids were optimized to achieve the best fitting of 
experimental data, i.e. the 2D position of the joints throughout the movements. The results showed that 
the models were able to reproduce accurately both movements. The best fitting was obtained with pointing 
tasks. The lower accuracy observed for squat-jumps may be explained by the transverse movement of joints 
which were not taken into account by the planar model. Furthermore, the pointing tasks would best meet the 
rigid body assumption compared to squat-jumps, for which the trunk was considered as a single limb. This 
study provides a new framework to model planar movements with a small number of meaningful kinematic 
parameters, allowing a continuous description of the kinematics and kinetics of the tasks. Further researches 
should investigate the implication of the control parameters in relation to motor control and validate this 
approach for three dimensional movements. 



1. Introduction 

Animal movement results from activation of complex musculo-skeletal systems whose architecture should 
be considered as the resolution, or part of resolution of environment's problems |Ber67j . Thus, they are very 
complex and potentially infinite. To cope with this complexity, human movements can be classified arbitrarily 
into four classes: 

- posture movements for which the subjects have to maintain a standing position; 

- pointing and reaching tasks for which the subjects have to move the distal extremity of a segment (i.e. 
the arm) to a target and eventually grasp it; 

- cyclic movements such as walking, running, cycling or swimming, for which the subjects have to move 
their body center of mass through cyclic behavior; 

- explosive movement for which the subjects have to move their body center of mass as fast as possible 
as in vertical jumping. 

In the literature, these movements have been analyzed through different paradigms and theories. However, 
even if they differ in control modalities by the central nervous system, they are, from a peripheral point of view, 
similar in the sense that they all imply joint displacements. When planning and executing such movements, 
the motor system must take into account several constraints in the gravitational field |IS89| . Considering any 
joint rotation, the more it approaches its maximal extension, the less the transformation of the segment angular 
velocity into linear velocity of its BMC is effective (geometrical constraint). Furthermore, it is necessary to 
decelerate the joint displacement before its maximal amplitude (anatomical constraint) in order to protect this 
joint from any damage. 
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Regarding to kinematics, the anatomical constraint implies that joint angular time history should match 
a sigmoid shape and thus a bell-shaped velocity profile |SL81| . which accounts for synergistic actuators' ac- 
tivations at a joint, i.e. agonist and antagonist muscle-tendon systems. The asymmetry of the joint angular 
velocity profile has been shown to depend of time to peak velocity and its intensity |ZSG86| . It was shown that 
the asymmetry increases together with the spatial accuracy requirements |MI90| and that faster movements 
are characterized by a longer acceleration phase |Pla95b| . 

Actually, one challenge in robotics and video games is to propose physic and behavioral simulations of 
human and animal subjects for the purpose of animating virtual characters or mechanical devices interacting 
with the environment. Formally, pointing an object in the environment requires the optimization of the inter- 
articular coordinations to move the distal extremity of the arm (i.e. the finger) throughout it's workspace 
[BLMlOl IBLM07] to a final posture which depends upon a complex combination of kinematic and dynamic 
factors [NMF99[ ISBHF95| . The first step in multi-joint movement simulation is to define a simple general 
heuristic model for motor control of mono- articular displacement. A movement can be considered as a transition 
from an initial to a final equilibrium configuration or posture |AAG"'"dT] . At an equilibrium configuration of 
the system, all forces and torques are balanced so that the system remain motionless. In response to deflection 
from the equilibrium configuration, the system generates additional forces tending to restore the equilibrium, 
if stability requirements are met. Theses authors concluded according to |GP71| that the equilibrium position 
of any physical system is not defined by output variables like EMG and forces but independently of them, by 
the system's parameters or control variables. 

During the recent years, the investigation of the processes implied in the control of movement planing 
and execution has made large use of simulation and optimization experiments. Using such approach, the 
present study aims at evaluating if movement execution can be controlled through a reduced set of kinematic 
meaningful parameters. According to the results provided by earlier researches, it was chosen to model the 
joints kinematics by using a flexible sigmoid shaped model of joint angle time history. For this purpose, 
pointing tasks and an explosive movement, the squat-jump, were modeled using three different sigmoid models 
based on kinematic parameters. Firstly, three kinematic based models of joint angle history meeting previously 
mentioned requirements are described. Then, these models were used to achieve best fitting of experimental 
data obtained from two specific movements, i.e. pointing tasks and squat-jumps. 

2. Methods 

2.1. Model of joint displacement 

2.1.1. General model. Accounting for a monotone evolution of a given angle and considering the anatomical 
constraint requirements, it is assumed that each angle 9 is characterized by the following properties (figure [T]): 

- at the beginning and at the end of the movement, the velocity and the acceleration are equal to zero; 

- the angle increases (respectively decreases) throughout the whole movement; 

- during the movement, the velocity increases (respectively decreases) until it reaches its maximum (re- 
spectively minimum), then decreases (respectively increases). 

More precisely, we try to determine a function from [0, T] to K of class . Let th,to, te be three instants 
such that 



<th <to <tc <T. 



(2.1) 



Let 0]^,9o, 9(. be three real numbers such that 



6'b < 6*0 < 9c or > 6*0 > 9^ 



(2.2) 



We assume that 

- 9 is constant and equals 9-^, on [0, <b] ; 

- 9 is constant and equals 9c on [tc, T] ; 

- there exists e G { — 1, 1} such that e9 is strictly increasing on [ib,^o] ; 



MODEL OF JOINT DISPLACEMENT USING SIGMOID FUNCTION. EXPERIMENTAL APPROACH 



angle 




time 



th to ^ te 
velocity 



K 













th to 


te 



time 



acceleration 











tb to \ 


(te 



time 



Figure 1. Shape of used sigmoid: angle, velocity and acceleration versus time (for the 
creasing case). 



- 60 is strictly convex on {tb,to) ; 

- e6 is strictly concave on (fo,ie)- 
We set 



Let K be the number defined by 



Since 9 is of class C^, we have 



£ = Sign {6, -Oh) G {-1,1}. 



K : 



max 9' it) if e = 1, 

te[tb,ic] 

min 9'{t) ife = -l. 
te[fb,te] 



e{th) = 9k, 9\th) = 0, 9"itk) = 0, 
e{t,) = 9,, 9\U) = 0, 9"{U) = 0, 
9{to)^9o, 9'{to) = K, 9"{to) = 0, 
yte{tb,to), e9"{t)>0, 
ytG{to,te), e0"{t)<O. 
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We consider a, /? € (0, 1) and k eR defined by 



a=*-^, /3 = §^, k = K^f^. (2.6) 

tc — ib — Oh Oc — Oij 

Applying the following change of scale, 

Vte[ib,ic], ii=^^e[0,l], (2.7a) 

To — lb 

yu e [0, 1], giu) = / > . (2.7b) 

the problem can be reformulated as follows: we look for a function g of class defined on [0, 1] satisfying: 

g(0) = 0, ,g'(0) = 0, g"(0) = 0, (2.8a) 

g(l) = l, 5'(1) = 0, 5"(1) = 0, (2.8b) 

g(a)=/3, g\a) = K = 0, (2.8c) 

VMe(0,a), 5"(t)>0, (2.8d) 

Vu e (a, 1), 5"(t) < 0. (2.8e) 

Remark 2.1. Under the assumptions given in (|2.8p we have necessarily (see |Cre09[ |BC12| ) 

k > max ^) > 1. (2.9) 



ift<ib, 



1 — a 

Finally, the function 9 is defined for all t £ [0, T], by 

9b, 



0(i)=<'(0e-eb)<?(r^^)+eb, iftb<t<ic, (2.10) 



— ^b 



if t > 



9c), 



This function is defined by 7 independent parameters: 

- 2 time scale parameters (tb and io 

- 2 angle scale parameters (0b and 

- and 3 shape parameters (a, /3, fc). 

Thus, 9 can be written under the form 9tt,,t^,9b,e^,a,i3,k- 

In the literature, there exist many sigmoidal functions. However, these models can not be used to solve 
because the sigmoids 

- are symmetric |MMMR96l IKumOOl IJNar97l IDG061 IYK03I IKS96] , 

- are defined on R and can not be used on a bounded interval |Dra95l ISC03I ILG63[ IBra70l IPea91| , 

- are characterized by not enough parameters (Fin47[ IFin52l IDeb79| . 

To our knowledge, there is not in the literature, a general construction of non symmetric sigmoid satisfying 
I, of class or C°° . 

In the field of movement analysis, R. Plamondon |Pla95al IPla95bl IPla98l IPCF03 j used a log-normal 
function 

^--<"^;^»-(-^^^^^^) 

to describe general movements and applied it to the Fitts task |MF541 IMFP64| . By using the Central Limit 
Theorem, he proved that, for a large number of agonist and antagonist muscles acting, the proposed function 
can model the behavior of the system. However, this work can not be applied to solve (j2.8p . The major concern 
with Plamondon's function remains in its asymptotic behavior at the end of the movement. Especially, it should 
be observed that the velocity tends to zero as t approaches +oo and then, the end of movement is not clearly 
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defined. Other works related to the log-normal law |LG63[ iBraTOI lPea91| do not solve this problem. Moreover, 
Plamandon's model contains not enough parameters to allow the solving of (|2.8p . 

The simplest idea to solve (|2.8p would be to apply Hermite's polynomial interpolation (e.g. [CdB81| ). but 
it can be showed that this method can not be used (see |BC12| 1. 

To allow the solving of the system (j2.8p . the model has to include three control parameters which have 
to be determined in relation to a, /3 and k. In the next section, three sigmoid models, SYM, NORM and 
INVEXP, previously described in |Cre09| . meeting the mentioned requirements are presented. These models 
were successfully used for both pointing tasks |Vil081 IVBML08I ILCB11| and squat-jump |Cre09l ICBL09| . 
Exhaustive theoretical description of the models will be given in a future paper [BC12] . 

2.1.2. The SYM model. 

The SYM model was built using a pseudo-symmetry approach. Its function g is defined by the three 
parameters a,/3 G (0, 1) and k > 1. 

Let ga.i3,k be a function of class from [0, a] to R satisfying (|2.8ap . (|2.8cp . and (j2.8dp . If the function g 
is defined from [0, 1] to R by, 

5(.) = l^-^''^^"^' (2.12) 
[l - 5ii_a4_/3^fe(l - u), if u > a, 

then, g is of class on [0, 1] and (|2.8p holds. Considering the function defined on [0, a] for all a,b > 

and K > 2 by 

F(«.f'^«)(u) = a (^1 , (2.13) 
a, b and k have to be determined so that (|2.8ap . (|2.8cp and (j2.8dp hold. We set 

For all (a,/?) G (0,1)^, for all k such that k > r^P/a, there exist (a, (2, oo) such that (|2.8al) . 

(p:Bc)) and ([2:8d| hold for function a, b and k stiU need to be defined. We set 



/3 

7 = 

It exists a unique X E (1/2, 1) such that 

X 



= — e (o,ei -l) . (2.15a) 
ka \ / 



e^-1) 1-^=7, (2.15b) 



and it follows 

1 /? , X 

K = , a = TT, — — . 

1-X' 1-e-^ a'^ 

By setting (a, &, k) ~ Gict, /3, fc), the function g is defined for all u G [0, 1] by 



(2.15c) 



9{u) = I . ' ^ " 2.16 

2.1.3. The NORM model. 

The NORM model (named from its relation to the normal law) function g is defined by three parameters 
a G (0, 1), p > and s > 0. 

We recall that the density function of the the normal (or Gaussian) distribution with mean m and variance 
is given by: 
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Considering the erf function defined by 

Vx e M, er: 



VT^ Jo 



e"* dt, 



the cumulative distribution function of the normal law is given by 



Vx e R, ^{x) = ^ erf 



X — m \ 1 
" 2' 



\/2s 



l-iiP 



For all p > 0, we define the bijection G from (0, 1) to R by 

Vue(0,l), G(u) = ln 
Finally, the function g is defined by 

Vie (0,1), g{t) = ^G{t)), 
<?(0) = 0, 
5(1) = 1. 

with a = G~^{m). 



(2.18) 



(2.19) 



(2.20) 



(2.21a) 
(2.21b) 
(2.21c) 



2.14. The INVEXP model. 

The INVEXP model (derived from the inverse exponential) function g is defined by three parameters 
A, > and a e M. For all a, for all A, ji, we set 

A 



a 

and we consider function ga^a defined by if a = 

9o.,ot ~ 

and if a > 



A + At 



e(o,i). 



(2.22) 
(2.23a) 



e [0, a), ga,a{y) = l - exp 



a{t — a) 



(2.23b) 



Vy e [a, 1], ga.aiy) = 1- 
For all a e M and for all a G (0, 1), we consider the function Ga,a defined by 



ifa>0, Ga,a=ga,a, 

ifa<0, Ga,a = g-a,l-a{^- ■)■ 



For all A, > 0, /a,u is defined by 



Vie (0,1), /A,;.(i)=exp 

/a,,.(0) = 0, 
/a,m(1) = 1- 
For all a G R A, /Lt > 0, /iA,/i,o is defined by 

hx,fi,a = f\,iJ.Ga,\/{\+ij,) 



(2.24) 

(2.25a) 

(2.25b) 
(2.25c) 

(2.26) 
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and finally the function g is defined by 



L 



{u)du 



Vie [0,1], g{t) 



L 



1 



(2.27) 



{u)du 



2.1.5. Definition domains of the sigmoid models. 

Each of the three functions is defined by three parameters. We will prove in jBC12| that for all fc > 1, 
there exist a part Sk of (0, 1)^ such that for all (a, (3) G Sk, there exist at least one sigmoid of kind g satisfying 
(|2.8p whose parameters can be determined by splitting (j2.8p in three non-linear equations which can be solved 
with a numerical solver. This part Sk is different for the three sigmoid models. The bigger part is obtained 
with the INVEXP model and is given by 



This domain is a polygonal part of [0, 1]^. The domain of a function satisfying (j2.8p can not be bigger thanks 
to (|2.9p . The domain Sf^^ of SYM sigmoid, which is also a polygonal part of [0, 1]^ is given by 



where rg is defined by (|2.14p . The domain 5'NORM NORM sigmoid can be determined numerically. It should 
be noticed that these three domains are symmetric according to the point (1/2, 1/2) and this point belongs to 
the three domains (Fig [5]). 

Examples of position and velocity curves obtained from the three models are provided in figure [31 Each 
sigmoid is defined by (ai,/3i), belonging to a fixed straight line and fc = 2. 

2.1.6. Specific properties of the sigmoid models. 

- The model SYM is very simple and fast to calculate; however, the class of this model is only versus 
C°° for the two other models. 

- Since the function erf is directly implemented in most numerical softwares, the model NORM is fast 
to calculate. However, the domain of this model has been obtained by symmetrization, and for some 
values of (a, /3), two different sigmoids can be obtained. Moreover, the determination of the definition 
domain is not trivial and the function does not meet the concavity and convexity requirements outside 
of it. 

- The part Sk of model INVEXP is the biggest, but this model is harder to calculate, because a numerical 
method of integration has to be used. Since efficient numerical methods exist, this problem can be 
overcomed and computation time remains reasonable. 

2.2. Experimental procedures 

2.2.1. Pointing task. 

22 subjects volunteered for this study. Nine right-handed male subjects (age, height) were asked to perform 
pointing tasks in the horizontal plane. Movements were performed for five directions and two distances (Fig 
U]). For each direction, two targets were placed at 60 and 80 cm from the shoulder. Directions of pointing task 
ranged regularly from 30° to 150° including pointing along the antero-posterior axis. The described position 
of the targets ensures that each of them is located inside the subjects workspace (Fig [S]) when considering a 
80 cm upper limb length and anthropometric data presented in |BLM10| . At the beginning of the movement, 
subjects had to position their arm so that the forefinger is located at 40 cm of the soulder in the antero- 
posterior direction. During the experiment, subjects sat on a chair whose height was adjusted so that the 
upper limb remained in the horizontal plane during the task and the trunk was immobilized by using straps. 
Video reflective markers were placed on the subjects at the shoulder (acromion), elbow (olecrane), wrist (middle 




(2.28) 




(2.29) 
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0.1 0,2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 OA 0.5 0.6 0.7 0.8 0.9 1 



(c) (d) 

Figure 2. Domain Sk of the three sigmoids for k — 1.2 (a), k = 1.5 (b), for fc = 3 (c), 
and for k = 4.5 (d). INVEXP, NORM and SYM domains are plotted in red, blue and green 
respectively. According to (|2.29p . SYM domain is empty for fc = 1.2 and k = 1.5. 



of radial and ulnar styloid processes) and forefinger extremity to allow further modeling of the upper limb. 
For each target, subjects performed three movements which were filmed at 25 Hz with a numeric camera JVC 
©Everio placed above the subjects and oriented vertically. Raw experimental data, i.e. the position of the 
joints throughout the movement, were extracted from videographic recordings. 

2. 2. 2. Squat jump. 

The squat-jump data was obtained from a previous work fB BMll| . Each of 13 other subjects performed 
10 vertical jumps. Instructions were given for keeping the hands on the hips during the movement to limit the 
contribution of the upper limbs to the performance. Furthermore, subjects were asked to do no countermove- 
ment. The jumps which did not meet both of these requirements were excluded from the study. In order to 
model the skeleton in a 4 rigid segments system, landmarks were placed on the left fifth metatarsophalangeal, 
lateral malleolus, lateral femoral epicondyle, greater trochanter and acromion. These landmarks define the 
foot, the shank, the thigh and the upper body (Head, Arms and Trunk: HAT). The subjects were filmed 
orthogonally to the sagittal plane at 100 Hz and the ground reaction force was recorded at 1000 Hz from an 
0R6- 7-2000 AMTI force plate. The center of mass (CoM) position of limbs was computed using anthropomet- 
ric data |Win90| . The whole body CoM position was determined on the one hand from kinematic data and on 
the other hand from force plate measurements using a double numerical integration procedure. For the latter, 
subject mass, initial body CoM position and velocity had to be set. These values were computed so that the 
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Figure 3. Examples of curves for angle and angular velocity, INVEXP model (a), NORM 
model (b), and SYM model (c). The boundaries of domain are plotted in dashed red line. 
Points {ai,f3i) are indicated by red circles. 



difference between CoM path obtained from kinetic and kinematic data was minimized in a least square sense. 
This optimization step was also used to synchronize both recording sources. 

2.3. Skeletal model 

For both tasks, the studied limbs were modeled as rigid bodies rotating around frictionless hinge joints. 
Given p limbs, the joint positions are defined by the points Aj (xj, yj) (j € {1.., . . . ,p}) with p = 3 and p — 4 
for pointing task and squat-jump respectively. From this definition, the position of the joints in the direct 
orthonormal reference frame (^0,i,j^ are related in the complex sense to limb lengths Ij and angles 6j (FigO 
and the affix of Aj is given by 

+ ^/„exp i^eJ (2.30) 

n=l \ k=l / 

for all j £ {2, . . . ,p}, where i is the imaginary unit and is the affix of ^i. 
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(a) : Upper view of the task environ- (b) : Targets (continuous line) and 
ment initial arm position (dashed line) 



Figure 4. Pointing task experience. 



0.2 - 



-0.4 - 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



Figure 5. Upper limb workspace |BLM10| 

It should be noticed that 9i and 9j{j € {2, . . . are segmental and joint angles respectively (FigurelH]). 
According to these definitions, the velocities and other derivatives of the joint positions with respect to time 
can be computed from corresponding derivatives of 9j and Ai . 

Remark 2.2. The determination of the positions, velocities and accelerations of the points A2, . . . , Ap require 
that the coordinates and further derivatives of Ai are known. Three cases should be considered to calculate 
xi and yi. 

- the point Ai is fixed; 

- the point Ai belongs to a simple curve such as a circle or a parabola; 

- the center of mass of the subject and angle j G 1, are known. 

2.4. Movement model 

Considering limbs as rigid bodies, the relative positions of joints are directly related to the angles 6j. The 
evolution of the angles throughout the movement was modeled using three types of sigmoid shaped curves. 
Theoretical results of this section have been already partially given in jCre09| and will be presented extensively 
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Figure 6. Skeletal model. 



in a future work |BC12| . We recall that to allow dynamic continuity, the required solutions should be models 
of class at least, defined on [0, 1] and satisfying (|2.8p . 

2.5. Data processing 

Experimental pointing task and jumping data were modeled using the sigmoid models. Specific procedures 
are described below. 

2.5.1. Shoulder path. 

0.22 
0.21 
0.2 
0.19 
0.18 
0.17 

0.86 0.87 0.88 0.89 0.9 0.91 0.92 0.93 

Figure 7. Shoulder trajectory: Experimental points and fitted circle. 
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Experimental data obtained from pointing tasks show that the shoulder joint path is well fitted by an arc 
of circle for all subjects (Fig. [7]). This movement was observed in spite of trunk immobilization. Since shoulder 
is a patella type articulation, this result is not surprising. The characteristics of the circle can be determined 
using a least squares method. The center f2 and radius R of the circle minimizing the distance sum of squares 
between experimental data {xk,yk) and theoretical data {Xk,Yk) were researched. Firstly, a direct method 
was used to minimize the sum of 

S = {{xk - Xkf + ivk - Ykf - R^f. 

Secondly, the sum 

S' - ^{xk - XkY + {vk - 

was minimized by using an iterative method. For this method, the results of the direct optimization were used as 
initial values for Q, and R. This final optimization was performed with the library Matlab least squares geometric 
element software, available at http : //www . eurometros . org/gen_report . p hp?category=distributionsfcpkey=14} 
For more details, see |Raz97[ IRaz98l lARWOli IEMM07[ IMK91I IFan90l ICjjg] . 

We set Aq = n, Iq = R and we consider then angle 6*0 defined by: 

00= e (-^,7r]. (2.31) 

We assume then that /q and Aq are constant and we obtain 

ZAi = ZAo + ^0 exp(i6'o). (2.32) 

where xo, yo and Zq are known. We add then Eq. (j2.30p : thus, assumption of remark l2. 21 holds. 

The mechanical system is plotted on Fig. |6l the coordinates of Aq and the lengths Iq, li, I2 and I3 are 
constant and the experimental data are then angle 6q, 61, 62 and ^3. 

2.5.2. Determination of sigmoid parameters. 

Sigmoid parameters were obtained from a multi-stage optimization procedure. First, t\^,tc,0\i,0c,Oi, P 
and K were estimated from experimental data. The scale parameters were defined so that the absolute peak 
angle velocity occurs between tb and <c and the angle velocity sign changes at the endpoints of this interval. 
Considering and 6^ as the angle values at tb and tc, the shape parameters a, /3 and k were determined thanks 
to Eq. (EH. 

First optimization consisted in minimizing the sum of square of differences between experimental angles 
9^ and those obtained from the sigmoid models cr* for each of n instants: 

n 

The optimization was achieved for each sigmoid with the Isqcurvef it function provided in Matlab software. 
Initial values of parameters were set from estimations of experimental data described previously. This opti- 
mization stage will be further referred to as local optimization. 

Secondly, differences between experimental and model reconstructed joint positions were minimized in a 
least square sense. The objective can thus be written as 

p n 

j=i i=i 

Compared to the previous stage, this optimization can be considered as global since for the latter, the 
parameters of the p sigmoids were determined simultaneously. Computation of model-based joint positions 
implies the lengths of the limbs to be provided. For the pointing tasks, the optimization was performed using 
(i) mean experimental limb lengths (semi-global optimization) and (ii) limb lengths as model parameters (global 
optimization). 
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2.5.3. Squat-jump specific procedure. 

The modeling of the jump focused on the position of the joints in a reference frame located at the distal 
extremity of the foot. Thus, the optimization consisted in fitting the experimental joint positions of ankle, 
knee, hip and shoulder with the model parameters in this reference frame. Since joints do not remain fully 
extended after the take-off, differences were not taken into account during the whole movement. This prevented 
the model from underestimating the necessary amplitude of joint extensions. Therefore, differences between 
experimental and model-based data were considered during the intervals corresponding to increase of vertical 
joint coordinates in the given reference frame (e.g. the error at the ankle joint was only taken into account 
while the vertical distance between the knee and the foot extremity increased). The objective S' was used to 
achieve this optimization stage, which will be further related as kinematic. 

Second stage of optimization included non-linear constraints on position, velocity and acceleration of the 
body CoM computed from sigmoid model. It was imposed that the body CoM position computed from both 
the sigmoid model and the force plate data were similar at the instant ti for which the marker located on 
the distal extremity of the foot started to move upward. At this instant, equality for the coordinates of both 
velocity and acceleration of body CoM obtained from kinetic and kinematic data was also required. Finally, 
body CoM vertical acceleration was constrained to be greater than -9.81 m.s^^ before ti ensuring that take-off 
occurs necessarily after ti. From ti to the end of the jump, the movement of was set so that kinetic 
and kinematic-based movement of the CoM were similar. This results in a continuous characterization of the 
movement position, velocity and acceleration. It should be noticed that using similar constraints for jerk and 
further derivatives could have led to description of class and higher. 

2.6. Modeling accuracy 

For each instant i and joint j, the optimization accuracy can be quantified by the difference between 
experimental data (xj and and sigmoid modeled data (Xj and Y^*): 

= ^{X]-x)f + {Y;-y]f. (2.33a) 

In further analysis, maximal Emax and mean £mean values of these differences are used to account for the fitting 
accuracy of the modeling procedures: 

Emax = maxe^.j, (2.33b) 
£mean = meaneij. (2.33c) 

2.7. Statistical analysis 

For both maximal and mean errors given in (j2.33bp and (j2.33c[) . the Shapiro- Wilk test reported unnormal 
distributions. Thus, statistical tests were realized on normally distributed logj^g of observations (i.e., errors and 
computation time). Firstly, anovas for repeated measures were performed for errors and computation time. 
When anovas reported significant results, post-hoc tests were performed to check for differences between the 
sigmoid models and the optimization procedures. All the tests were realized with <S [R Dll) and statistical 
significance was set at 95% confidence level, i.e. p < 0.05. 

3. Results 

The results obtained from the 304 pointing tasks and 120 squat-jumps are summarized in the tables [5] to 
[7] given in the appendix |21 Examples of experimental and modeled data are provided for the pointing tasks 
and squat-jumps in the figures 151 to [TT] and 1 141 to [T71 respectivelv. 

All tables are given in appendix [Xj 
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0.2 0.3 0.4 0.5 0. 



(a) : Angle 




0.1 0.2 0.3 OA 0.5 0.6 0.7 



(b) : Angle 1 




(c) : Angle 2 



(d) : Angle 3 



Figure 8. Angles (in degrees) according time. Experimental data 6*^ are plotted by black 
points, INVEXP sigmoid model is plotted in red continuous line, NORM sigmoid model is 
plotted in blue continuous line, SYM sigmoid model is plotted in green continuous line. On 
the figure, each points of coordinates (^t'^\0]^^^ and (^ti''\0i''^^ are plotted by a circle. We 

add also the tangent of the curves in 



^0 



3.1. Pointing tasW 

Basic descriptive statistics of measured values are given in tables [2] and [3l 

Remark that, if Q — 304 is the number of measures of error Si, we have by setting rji = logj^Q (e^) 

/ \ i/Q 

n^O = io(*^-"') = io'°sio(^), 

which means that the geometric mean of error is equal to 10™ where m is the mean of log^^g of error: see table 

SI 

Recall that for p G [0,1] 

- '***' means p < 0.001; 

- means p < 0.01; 

- means p < 0.05; 

- '.' means p < .1. 
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Figure 9. Angles (in degrees) and angular velocity, acceleration and jerk, for optimization 
for each angle; Experimental data 6** are plotted by black points, Smoothing data 9j are plotted 
by dashed black line, INVEXP sigmoid model is plotted in red continuous line, NORM sigmoid 
model is plotted in blue continuous line, SYM sigmoid model is plotted in green continuous 
line. 



Among the three anovas performed for computation time, maximal and mean errors, the highest p- value 
was equal to 1.48e — 119 (***) suggesting that both optimization methods and sigmoid models are associated 
to significantly difTerent results. 

Considering both mean and maximal errors, post-hoc tests revealed greater adequation of original data 
with semi-global optimization procedure than with local one {p = 0) for each of the three sigmoid models. 

Computation time reported for semi-global optimization was significantly higher than durations obtained 
with local method for SYM and INVEXP models (p = 0) but not for NORM one (p = 1). 

Comparing semi-global and global methods revealed no significant difference for computation time {p — 1), 
maximal {p > 0.3291) and mean {p > 0.6166) errors. 

For the local optimization case, no significant difference was found between the three sigmoid models 
when comparing maximal {p > 0.7999) and mean {p > 0.7737) errors. Similar results are obtained for global 
(maximal: p > 0.1314, mean: p > 0.1957) and semi-global (maximal: p > 0.1569, mean: p > 0.1878) cases. 
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Figure 10. Angles (in degrees) and angular velocity, acceleration and jerk, for global opti- 
mization; Experimental data 6j are plotted by black points. Smoothing data 9j are plotted by 
dashed black line, INVEXP sigmoid model is plotted in red continuous line, NORM sigmoid 
model is plotted in blue continuous line, SYM sigmoid model is plotted in green continuous 
line. 



3.2. Squat Jump 

Basic descriptive statistics of measured values are given in tables [S] and [51 Examples of experimental and 
modeled joint angles time histories are presented for both kinematic and dynamic optimization methods in the 
figures [HI and [T5] respectively Time histories of relative joints position are presented in the figure [TBI 

Among the three anovas performed for computation time, maximal and mean errors, the highest p- value 
was equal to 1.856e— 136 (***) suggesting that both optimization methods and sigmoid models are associated 
to significantly different results. 

Compared to constrained optimization, unconstrained method executed faster {p — 0) and fitted better 
kinetic original data for both maximal and mean errors (p = 0) whatever the sigmoid model. 

Considering the unconstrained optimization, no significant difference was found between the sigmoid mod- 
els for maximal (j) > 0.06547) and mean {p > 0.06293) errors. 

Post-hoc tests revealed significant differences between the sigmoid models for the constrained optimization 
method. Maximal errors mesaured for NORM model were higher than those obtained from SYM {p — 3.73e— 06 
(***)) and INVEXP (p = 0.02128 (*)). Considering mean errors, SYM model fitted best the original data 
compared to NORM {p = 3.282e - 10 (***)) and INVEXP (p 5.018e - 06 (***)). 
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Figure 11. Trajectories of points Ai, Ai, A3 and A^,. INVEXP sigmoid model is plotted in 
red continuous line (for global optimization) and red dashed line (for optimization for each 
angle). NORM sigmoid model is plotted in blue continuous line (for global optimization) and 
blue dashed line (for optimization for each angle). SYM sigmoid model is plotted in green 
continuous line (for global optimization) and green dashed line (for optimization for each 
angle) . 



4. Discussion 

This study evaluated different optimization methods to fit joint trajectories produced during pointing 
tasks and squat jumps. The evolution of joint angles during the movements was modeled using three sigmoid 
shaped functions. Assuming a constant length of the limbs, the whole movement was reconstructed from the 
sigmoid models parameters. For each movement type (i.e. pointing tasks and squat jumps) and sigmoid model, 
different optimization methods were investigated. 

In the literature, only Plamondon used a similar approach. However among the published articles, exper- 
imental data were presented only in |Pla98| . Furthermore, no quantitative results were provided and the data 
was presented for a single subject. This does not allow to compare the present models with Plamondon's one. 
However, as mentioned earlier, the models used in the present study are defined on a bounded time interval 
contrarily to the log-normal models for which the end of the movement is not clearly defined. 

4.1. Rigid bodies assumption 

Differences between original and reconstructed data were lower for pointing tasks than for squat-jumps. 
This result could be explained by the relatively greater amplitude of the joint trajectories during the jumping 
movement. Moreover, the modeling of the skeleton assumes rigid bodies between the joints. Considering the 
pointing tasks, it can be supposed that the length of the modeled limbs is quite constant. This assumption 
is supported by the similarity of the errors observed for global and semi-global methods. The rigid bodies 
assumption would be less true for squat-jump, especially for the trunk limb. Indeed, the spine is composed 
of many joints which allow bending of the trunk and thus, the trunk may be divided into two |KdLT^96[ 
ldLKBT92l IPGD96| or three |dL93| segments to ensure that the rigid bodies model is close enough to the 
reality of the movement. 
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Figure 12. Histograms of Smax for INVEXP, NORM and SYM sigmoid models, for semi- 
global optimization Continuous and dashed lines indicate mean and 95 % quantile respectively. 



4.2. Planar movement 



The modeling methods proposed in this study deal with planar movements. The higher errors obtained 
with modeling of squat jumps may be explained by the movement of the joints along the transverse axis, 
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Figure 13. Histograms of logiQ(£max) for INVEXP, NORM and SYM sigmoid models, for 
semi-global optimization Continuous and dashed lines indicate mean and 95 % quantile re- 
spectively. 



especially for the knee. In comparison, pointing tasks would be closer to a real planar movement since the 
movement is performed on a planar surface. 
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Figure 14. Time histories of joint angles. Dotted and plain curves correspond to experimental 
and kinematic stage modeled data respectively. Vertical lines indicate ti and take-off instants. 



4.3. Optimization methods computing velocity 

Considering the pointing tasks, computation lasted longer for semi-global method than for local one. 
Global optimization executed with similar velocity compared to semi-global method. Thus, global optimization 
should be used unless specific purposes are researched. For the squat-jumps, the present results show that 
unsurprisingly, using the constrained method is much more longer than the unconstrained optimization. 

4.4. Optimization methods accuracy 

For pointing tasks, accuracy of the model was higher for semi-global optimization than for local one. 
This suggest that modeling should consider the joints movements together to achieve better fitting of original 
data. It should be noticed that local optimization could have considered the dependence of the distal joints 
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Figure 15. Timehistoriesof joint angles. Dotted and plain curves correspond to experimental 
and dynamic stage modeled data respectively. Vertical lines indicate ti and take-off instants. 



trajectories to the proximal ones. Global optimization did not lead to better results compared to semi-global 
method. This result is consistent with both the planarity of movement and rigid bodies assumptions. 

4.5. Sigmoid models accuracy 

For both pointing tasks and squat jumps, similar accuracy was obtained with the three models of sigmoids. 
Among the two movements and the optimization methods, it appears that the NORM model allows fastest 
computation. Considering SYM and INVEXP models, their relatively slower execution can be explained by 
the non-linear equation solving and the numerical integration respectively. NORM model formulation takes 
advantage of the native implementation of the erf function in Matlab software thus ensuring fast computation. 
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Figure 16. Time histories of joints relative positions. Dotted and plain curves correspond 
to experimental and kinematic stage modeled data respectively. Vertical lines indicate ti and 
take-off instants. 



4.6. Practical considerations 

The present results show that joint trajectories during planar movements such as pointing tasks or squat- 
jumps can be modeled using meaningful kinematic parameters. Table [T] presents a summary of the results 
obtained with the different optimization methods and sigmoid models. 
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Figure 17. Time histories of joints relative positions. Dotted and plain curves correspond 
to experimental and dynamic stage modeled data respectively. Vertical lines indicate ti and 
take-off instants. 

Among the three sigmoid models tested in this study, it appears that the NORM model is computed 
faster and allows better data fitting of the pointing tasks than other models. On the contrary, for squat-jumps, 
INVEXP and SYM models fitted better original data. From these results, it can be suggested that INVEXP 
and NORM models should be used preferentially. Indeed, the INVEXP model did not lead to better results 
and needs substantial computation time compared to other models. Despite the important computation time. 
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SYM 


NORM 


INVEXP 


Computation time 


+ 


+++ 


+ 


Definition space 


+- 


+ 


+++ 


Accuracy 











Mathematical regularity 


+ 


+++ 


+++ 



Table 1. Summary of the results obtained for the sigmoid models. + and indicate the 
existence or absence of advantages respectively 

INVEXP model may be useful for modeling specific movements, especially fast movements, which may not 
allow a good fitting with NORM model due to the relatively small definition domain of this model in the 
a, (3, K space. For relatively slow and smooth movements, NORM model should be primarily used. 

5. Conclusion 

This study shows that complex planar movements can be modeled accurately by using a small set of mean- 
ingful kinematic parameters defining the time history of joint angles. This approach can provide a continuous 
description of the movement and thus may be used to analyze the evolution throughout the movement of param- 
eters which need differentiation of raw data with respect to time without performing numerical computations. 
Especially, this could avoid well known magnification of error resulting from such procedure. Moreover, the 
flexibility of the new sigmoid models should lead to increased realism of movements obtained from procedural 
animation. Further researches should assess the relevance of such modeling strategy for three dimensional 
movements and the relation between the model parameters and the central nervous system processes implied 
in motor control. 
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Appendix A. Set of tables of statistical results 
A.l. Pointing task 



data 


method \ model 


SYM 


NORM 


INVEXP 


computation time 


local 


0.77() ± 0.24 


0.471 ± 0.23 


1.41 ± 0.17 


semi-global 


1.361 ±0.33 


0.309 ±0.32 


1.836 ±0.37 


global 


1.248 ±0.45 


0.079 ±0.41 


1.66 ±0.5 


maximal error 


local 


-1.915 ± 0.22 


-1.917±0.22 


-1.908±0.21 


semi-global 


-2.102 ±0.2 


-2.121 ±0.2 


-2.093 ±0.2 


global 


-2.118 ±0.2 


-2.139 ±0.2 


-2.109 ±0.2 


mean error 


local 


-2.415 ±0.22 


-2.422 ±0.21 


-2.412±0.21 


semi-global 


-2.591 ±0.18 


-2.612±0.18 


-2.586 ±0.18 


global 


-2.6 ±0.18 


-2.62 ±0.18 


-2.594 ±0.2 



Table 2. Elementary statistics on log^Q (mean ± standard deviation). 



data 


method \ model 


SYM 


NORM 


INVEXP 


maximal error 


local 


2.768 


2.894 


2.781 


semi-global 


1.86 


1.609 


1.725 


global 


1.645 


1.486 


1.602 


mean error 


local 


0.824 


0.824 


0.812 


semi-global 


0.521 


0.451 


0.484 


global 


0.493 


0.447 


0.47 



Table 3. 95 % cases error in centimeter. 



data 


method \ model 


SYM 


NORM 


INVEXP 


maximal error 


local 


1.2169 


1.2169 


1.237 


semi-global 


0.7912 


0.7567 


0.8076 


global 


0.7618 


0.7265 


0.7777 


mean error 


local 


0.3846 


0.3846 


0.3873 


semi-global 


0.2565 


0.2444 


0.2592 


global 


0.2513 


0.2401 


0.2545 



Table 4. Geometric mean of errors in centimeter. 
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A. 2. Squat Jumps 



data 


method \ model 


SYM 


NORM 


INVEXP 


computation time 


kinematic 


1.334 ±0.27 


0.26 ±0.24 


0.978 ±0.19 


dynamic 


2.114 ±0.28 


1.311 ±0.36 


2.096 ±0.24 


maximal error 


kinematic 


-1.311 ± 0.16 


-1.33 ± 0.17 


-1.297±0.16 


dynamic 


-1.033 ± 0.21 


-0.962 ± 0.2 


-1.001± 0.19 


mean error 


kinematic 


-1.771 ±0.13 


-1.784 ±0.14 


-1.76 ±0.13 


dynamic 


-1.606 ±0.16 


-1.535 ±0.18 


-1.554 ±0.14 



Table 5. Elementary statistics on log^g (mean ± standard deviation). 



data 


method \ model 


SYM 


NORM 


INVEXP 


maximal error 


kinematic 


8.392 


8.492 


8.492 




dynamic 


21.134 


26.257 


21.213 


mean error 


kinematic 


2.688 


2.688 


2.882 




dynamic 


4.935 


5.842 


5.034 



Table 6. 95 % cases error in centimeter. 



data 


method \ model 


SYM 


NORM 


INVEXP 


maximal error 


kinematic 


4.8902 


4.8902 


5.047 


dynamic 


9.2603 


10.9063 


9.9664 


mean error 


kinematic 


1.6946 


1.6946 


1.738 


dynamic 


2.4802 


2.9156 


2.7936 



Table 7. Geometric mean of errors in centimeter. 
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